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Abstract 

We recall the origin of angular ordering of soft parton emission in the region of 
small X and show that this coherent structure can be detected in associated dis- 
tributions. For structure functions at small x and at fixed transverse momentum 
the angular ordering is masked because of the complete inclusive cancellations of 
coUinear singularities for a; — *• 0. Therefore, in this case the dependence on the 
hard scale is lost and the angular ordered region becomes equivalent to multi- 
Regge region in which all transverse momenta are of the same order. In this 
limit one derives the BFKL equation. In general such a complete cancellation 
does not hold for the associated distributions at small x. The calculation, which 
requires an analysis without any collinear approximation, is done by extending 
to small X the soft gluon factorization techniques largely uses in the region of 
large x. Since one finds angular ordering in the both regions of small and large 
X, one can formulate a unified evolution equation for the structure function, a 
unified coherent branching and jet algorithm which allows the calculation of as- 
sociated distributions in all x regions. Such a unified formulation valid for all x is 
presented and compared with usual treatments. In particular we show that the 
associated distributions at small x are sensitive to coherence. By replacing angu- 
lar ordering with the multi-Regge region one neglects large singular contributions 
in the associated distributions. 
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1 Introduction 



Soft gluons emission in perturbative QCD takes place into angular ordered regions [IP-IQ. 
Since this coherent emission factorizes, one can extend the branching process, which sums 
collinear singularities to leading (and next-to-leading) order, to include angular ordering. 
This coherent branching process allows one to compute parton distributions either by nu- 
merical methods via Monte Carlo simulations or by analytical methods via "jet calculus" 
algorithms [|l], Experiments are well compatible with coherence and for a recent review 
see Ref. 0. An important case in which soft gluons are involved is deep inelastic scattering 
(DIS) in the region of large x with x = Q'^/2pq and = —q^ (see fig. 1). Coherence 
in the structure function F{x, Q) for x —>■ 1 can be taken into account in the Altarelli, 
Parisi, Gribov, Lipatov, Dokshitzer (AP) equation |^ by using as evolution variable the 
angle of emitted partons. Coherent branching at large x reproduces |§] the leading and 
next-to-leading contribution of the anomalous dimension 7a; (a^) 

l.ias) = Q' ^\^^}^\ FUQ) = £dxx-F{x,Q), (1) 

in the limit — oo which corresponds to the region of large x. 

The small x region in DIS appears quite different from the one of large x (or finite x in 
general). For x — > one focus the attention on the soft exchanged gluons ki, ■ ■ ■ ,kr rather 
than on the emitted ones pi,---,Pr- The basis for the analysis of the small x structure 
function is the Balitskii, Fadin, Kuraev and Lipatov (BFKL) equation ^ which was found 
almost two decades ago in the field of soft physics to discuss the rise of the total cross section. 
The dominant phase space region considered is the so called " multi-Regge" region in which 
subenergies are large and all emitted transverse momenta are of the same order of magni- 
tude. In general, if the ratio of transverse momenta vanishes, one has collinear singularities. 
However for fully inclusive distributions all these collinear singularities completely cancel, so 
that the regions of very different transverse momenta are not particularly important. 



Later it was realized [^, [T^ that this equation could be extended to describe the structure 
function for x 0. This equation involves the unintegrated structure function J-'{x, k^) 
related to the structure function by 

F^{Q) = j d^kT^{k')Q{Q-k), k = \\i\, (2) 

where k is the total transverse momentum of emitted gluons, k,. = I]Pfj. For ^ 
one probes the small x region and in this limit the k-structure function satisfies the BFKL 
equation 



J-.(k) = 6\\^) + ^ ffPl [j-^(k + pj - e{k - pt) ^(k)] , as ^ ^ , (3) 

UJ J npt TT 



with Ca_ = 3. The corresponding — > anomalous dimension is given by an implicit 
equation and has an expansion in a^g/uj^. The first six terms are 

iMs) = ^ + 2C(3) ^ + 2C(5) § + ■■■, (4) 
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with ({3) = 1.202 . . . and ({5) = 1.037. . . the Riemann functions. For the expansion up to 
14 terms see |0. As well known 7(j(as) has a singularity at a; = u;* = 4 In 20:5 which implies 
that the structure function rises asymptotically at small x. The two terms in the kernel of 
(0) have a simple interpretation. The first is positive and is the contribution for the emission 
of a gluon of transverse momentum p^. Thus, before the emission, the k-structure function 
has momentum k + Pj. The second negative contribution in the kernel corresponds to the 
virtual correction in which the virtual transverse momenta p^ is bounded by the total emitted 
transverse momentum k. There are two related features which make the BFKL equation 
apparently very different from the AP equation. The first is that there is no any ordering 
in the real emission and virtual momenta. The second is that the collinear singularities for 
Pt — > present in the real and virtual contributions cancel completely. This last feature can 
be further understood if we consider the exclusive representation of the k-structure function 
in eq. (^. To this end we need to introduce a cutoff fi both for the real and the virtual 
contribution and the BFKL equation becomes equivalent to the exclusive representation (see 
the kinematical diagram in Fig. 2) 

^^(k) = 5'\k) +Y. L li^ dz,z^^ Ar{z,, h) 5^(k - k,) , k, = p„ + k,_i 




Ai?(zi, ki) = exp <j -as I — / ' J> = exp -as In — In % e{kti - , k 




(5) 

The function ^^{zi, ki) sums all virtual contributions and corresponds to the so called gluon 
"Regge trajectory". In this exclusive representation we have included, together with the 
collinear cutoff /i, also the upper limit for the emitted transverse momenta given by the hard 
scale Q (see Due to the complete cancellation of collinear singularities at the fully 

inclusive level and to the scale invariance of the measure (Ppt/pf, the /i and Q dependence in 
the k-structure function (|^) enters only as higher twist corrections, i.e. as powers of fi^/Q'^ 
and for /x^/Q^ the representation (||) is equivalent to the BFKL equation (^. Since 
Q enters only to higher twists one has that for x — the hard scale is lost to the present 
level of accuracy. In DIS for x ^ the hard scale Q is recovered only through the phase 
space boundary in the k dependence of the coefficient function via the high energy (or k) 
factorization (see |TT|). On the other hand to the present accuracy the hard scale for the 
argument of as can not be identified. Since only the coefficients of a^jiJ^ are known, a scale 
change between Q and Q' 

asm ^ MQ) ^ «|(Q) 26 In(^) + . . . , 6 = (IIC^ - 2n,)/12vr 
iij uj u) Q' 

can not be discriminated. The control of the hard scale can be obtained by the calculation of 
subleading corrections as (Jsjoj^ in the anomalous dimension. The first steps in this direction 



are reported in Ref. |12 



From the above discussion it seems that the dynamical features of the two DIS regions of 
small and large x is quite different . However it was shown that the two regions have the same 
structure of coherence. The calculation of the structure function for x ^ was done in 
the framework of hard processes thus considering not only the soft singularities for x ^ but 
also the collinear singularities when the ration of two emitted transverse momenta vanishes. 
This region is beyond the multi-Regge region. The result was the discovery of coherence 
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for x — >■ 0. One finds that for x — initial state soft gluons are emitted into the angular 
ordered region as in the case of large x. In this way the description of DIS is unique for all 
region of x. 

Notice the accuracy required for this calculation. Since in fully inclusive distributions all 
coUinear singularities cancel for x — > 0, one needs to perform a calculation without making 
any coUinear approximation. In the light-cone expansion infinite leading twists operators 
contribute and one has to make a calculation to all loop accuracy. Actually, by using the 
soft gluon factorization techniques for x 0, it was shown [^] that it is possible to attain 
such an accuracy by appropriate use of soft gluon coherence and of cancellation of soft 
singularities between real-virtual contributions. One can use soft gluon techniques also for 
X — because this region is dominated by soft exchanged gluons ki, - ■ ■ ,kr and this implies 
that, a part for the first fast gluon pi, all remaining emitted gluons are soft. 

Many data are now available from Hera in the small x region and then the analysis of 
coherence at small x can be attempted. To help this analysis, in this paper we present 
a systematic account of coherence in all x regions. In particular, by using the coherent 
branching process discussed above we show how computes the associated distributions via 
a generalization of jet calculus This allows one to study quantitatively the role of 

coherence for x — and to single out the differences between the present formulation with 
coherence and the exclusive representation in (H) or the usual AP evolution equation naively 
extended to small x. 

The result of this analysis in the case of small x 0-0] can be summarized in terms of 
very simple extension of the exclusive representation of the BFKL k-structure function in 
It involves the following two modifications: 

i) the emission takes place in the angular ordered region. When all exchanged gluons are 
soft angular ordering can be expressed by including a factor Q{pti — Zi-ipti-i) in the phase 
space for each emission in (|^). The coUinear cutoff is needed only for the first emission with 
z = 1, i.e. pti > ZoPto = (this cutoff is not needed for the k-structure function); 

ii) also the virtual corrections involve angular ordering. The form factors A/j(zj, kj) in (^) 
are modified by substituting the cutoff pt > fi hj the angular ordered constraint pt > zpu- 

On the basis of this branching structure one deduces [0, Q an evolution equation for the 
k-structure function at small x which involves as evolution an angular variable. From this 
equation one observes that there is a complete cancellation of the coUinear singularities. In 
this case angular ordering becomes equivalent to the multi-Regge region (^. On the contrary 
we show that in the associated distributions, such as the associated multiplicity, the coUinear 
singularities do not cancel. In this case one can not use the multi-Regge region in (|) but 
must take into account angular ordering and the above modifications. 

The fact that both for small and large x parton emission and virtual corrections factorizes 
and have the same structure of coherence, i.e. angular ordering, allows one to deduce a 
unified branching process for all x regions. This has been used to construct a Monte Carlo 
program ITBI including coherence in the small x region which has been used to study the 



differences with the one loop formulation and to compute numerically DIS and heavy flavour 



production processes. Moreover one can write |T^ a unified evolution equation for the k- 
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structure function which extends the AP equation to include the coherence at small x 0, 
We show that this equation, which x —>■ reduces to the BFKL equation, can be used as the 
bases for the jet calculus to compute the associated distributions. 

In Sec. 2 we recall the results of Refs. of the DIS analysis for x — > together 

with a; ^ 1. We describe the accuracy of the approximations in the various region of x. In 
Sect. 4 we discuss the branching process with coherence for all x and in Sect. 5 we recall the 



unified evolution equation |jTj] for the k-structure function valid in all x regions. In Se. 6 
we discuss the small x region. At the fully inclusive level we show that from the evolution 
equation with coherence one deduces the BFKL equation valid in the multi-Regge region. At 
the semi-inclusive level (associated distributions) we show that the present formulation with 
coherence can not be approximated by using only the multi-Regge region as in the BFKL 
equation. In Sect. 7 we describe the extension of jet calculus to include the coherence at 
small X and compare with a formulation in the multi-Regge region. Finally Sect. 8 contains 
some conclusion. 



2 Derivation of recurrence relation and solution 



We consider deep inelastic scattering represented at parton level in Fig. 1 where q is the 
hard colour singlet probe and p' represents the recoiling system of partons, such as, for 
instance, a pair of heavy quark-antiquark. For simplicity we assume that the incoming p 
and the outgoing partons pi...pr are gluons. We take p = {E,0,0,E) and introduce a 
light-like vector p = {E, 0,0 — E) in the opposite z— direction. For the emission momenta 
we use the Sudakov parametrization 

— 2 

Pi = ViP + ViP + Pu, ii = — = , s = 2pp, (6) 

where In^j is the (pseudo-)rapidity. For the exchanged momentum one has k = xp + xp + kt. 
We assume the frame q = p — xp with x = Q"^ /2pq and consider = — large. In this 
limit we are interested in the region where all yi are small. 

The cross section of this process factorizes into the elementary hard scattering distribution 



and the structure function. For x ^ this factorization can be generalized [|T^ to include 
higher order corrections by introducing the structure function at fixed k. To leading twist, 
in the gauge ri = p = q + xp^p', the amplitude of Fig. 1 can be factorized Q as 
in Fig. 3 into the elementary vertex V'^(g, k) and the multi-gluon emission subamplitude 
m|[)(p, k;pi, . . .pr) which is a matrix in the colour indices of gluons p,pi, ■ ■ - Pr and k. In 
this gauge, for x ^ 0, the elementary vertex V^, is proportional to and one defines the 
multi-gluon subamplitudes M.r(pi, . . . = (p'^/2pp) m.^^\p, k;pi, . . . ,pr). One introduces 
the structure function at fixed total transverse momentum 

k^ k) = - 5] i / flidpi) \Mr{pi... Pr) I' 5(1 - X - F,)(52(k - k,) , 

d^p TV \n 
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with I Mr{pi . . .pr) f the spin and colour average of the multi-gluon squared subamphtude. 
The structure function is obtained by integrating over k 

F{x,Q) = I d^ktJ^{x,k,Q), T{x,k,Q) = QiQ-k)F{x,k), A; = | k | . (8) 

By integrating over i/i in (J^) one finds infrared (IR) divergences which are cancelled by 
virtual contributions. These cancellations are most easily seen in old fashion perturbation 
theory in which virtual gluons are on-shell. To order a^, the distribution | Mr \ is given 
by integrals over i = n — r on-shell virtual momenta Vi - ■ ■ ve. We introduce the on-shell 
momenta {gi, . . . , g„} which includes both the r real {pi . . . ,Pr} and the i = n — r on- 
shell virtual {f i . . . v^} momenta and we use the Sudakov parametrization = yip + yip + 
qti with yi = qu/yiS. One has to consider contributions with various ordering in the yi 
variables of the on-shell real and virtual momenta and for a given contribution we denote 
by fr,e {pi ■ ■ -Pr'jVi ■ ■ - Ve) the corresponding integrand in which the spin and colour sum is 
already performed. We then introduce the index i?„ which identifies the real {pi} or virtual 
{vi} momenta in the set {qi} and denote by fn^iqi ■ ■ - qn) the corresponding integrand. We 
then have then the following representation of the structure function. 

oo „ n 

X^(x,k,Q) = 5:^?|"/nW0(Q-fe)0L..n E/«n(gi---gn) 

n=l i=l (9) 

y,S{l-x-YnJ 6' (k - knj , Ol,^ ^ d{y, - y,) ■ ■ ■ - y„) , 

where Yr^ = and k^^^ = k^ are the contributions from the real gluons and are defined 
in (1^). Obviously the integrand . . . ^n) depends on the specific index R^. The bound 

Q > qu is coming from > k^ in (|^) and the ordering of the y^ variables is coming both 
from the statistical factor 1/r! for real momenta and from all possible ordering of the virtual 
on-shell momenta. 

We need now to compute ■ ■ - qn)- We will do this in the strongly y-ordered region 

1/n < • ■ • < 2/1 , (10) 

in which we can apply the soft gluon techniques and factorize the contribution of the softest 
gluon, either real or virtual. This region gives the leading contribution both for x — and 
for a; — >■ 1. The fact that for x — > 1 the strongly ordered phase space gives the dominant 
contribution is well known The fact that the ordering (|10D gives also the dominant 
contribution for x — is due to the infrared regularities of the inclusive distributions. Since 
in (^) the various singular contributions for ?/j — cancel, the integration regions yi < x give 
a contribution to the structure function of order x which can be neglected for x — 0. On the 
other hand in the remaining regions yi > x the real-virtual cancellations are inhibited by the 
presence of 5(1 — x — Yn^) in the integrand. Therefore x plays the role of an infrared cutoff 
and logarithmic contributions in x are generated by the lower bound in the yi integrations. 

In the remaining part of this section we shall recall ^, Q that by factorizing the softest 
(real or virtual) gluon g„ one deduces the following recurrence relation connecting f^^ and 
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/ Rn-1 

H fRniqi ■■■qn) gl (dq^) 5(1 - x - Yn,J6^(k - k^J ~ ^ /H„_i(gi • • • qn-i) "5 — — f 

X {<^(1 Yr„_, - Vn) 5'(k - k^_, + q„) 

- 5(1 - X - r^„_J 52(k - k^„_J [1 - 0(y„ - x)e(g„ - fc)] } . 

(11) 

As we shall recall, this recurrence relation is valid in the infrared limit (i.e. j/„ <^ Ui) for 
both X — > (without collinear approximation) and x — > 1 (in the collinear approximation). 
The first contribution in the curly bracket corresponds to the case in which the softest 
gluon is emitted, while the second corresponds to the case in which g„ is a virtual gluon. 
Consequently, the momentum of the softest gluon enters in the the 5— function only 
for the first contribution. This recurrence equation gives rise to an integral equation for 
the A;— structure function. One needs to fix the total transverse momentum k because this 
variable enters in the coefficient of the 5— function of the last term of ([TT|). 



2.1 Recurrence relation for 

We summarize here how to obtain from the soft gluon factorization the recurrence relation 
(|ll|) and discuss the relative accuracy. First we recall the case of large x and then we discuss 
the case of x both small and large. 

Large x. As well known the factorization of the softest (real or virtual) gluon g„ in this 
region gives the colour matrix factor 

jR„_i(gn)l^ I (5(1 - X - yR„_, - Vn) 5^(k - k^,„_, + q„) 

^ . (12) 

-5(l-x-F^_J52(k-k^„_j|. 

This result is valid in the soft approximation (?/„ <^ yi) and no collinear approximations 
are needed. In the first contribution g„ is the real emitted gluon thus it contributes to the 
5-functions in (H). In the second contribution the gluon g„ is virtual and does not contribute 
to the 5-functions. The soft emission factor is the squared of the eikonal current 

jR„_,(gn) = -— Tp+ ^T,^ + ^T,,, T,, = T,- Y T,., (13) 

with Tq. and Tp the colour matrices of external gluons, and Tp/ the colour matrix of 
the recoiling system p'. The emission factor is universal since it depends only on the fastest 
charges, therefore it does not depends on particular Feynman graph contributions. Moreover 
the factor is the same for g„ real and virtual. 

The current J is a matrix in the space of {p, qi, - ■ ■ , qn-i,p'} colour indices and in the use 
of (p!3| ) the difficulty is the colour algebra. Since the integrand is given by the average 
over the colour charges, the factorized expression in (12) gives a recurrence relation for f^^ 
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only if the colour factor in (12) is approximately proportional to the unit matrix. In this 
case of large x this is true only in the coUinear limit in which the emitted gluons qi are 
taken parallel to the incoming one, i.e. qi ~ yip. By using colour conservation in (0) and 



Tp = Ca = Nc, in this limit one has 

rH._.wr--c.f^-^)^^%. (14) 

J \pqn p'qnj q.t 

For large x the coherent emission is obtained in this approximation to the leading coUinear 
accuracy. Therefore the analysis in this case is to the leading IR and coUinear accuracy. 

Small and large x. This analysis can be extended to include the region of small x. Here 
the analysis seems much more difficult since one can not make any coUinear approximation. 
Moreover soft gluon factorization techniques are typically used in the case in which the softest 
gluon qn is softer than all other external q^ and internal lines. In the a; — ^ case instead we 
have that the exchanged gluon fc„ could be softer than g„. ActuaUy the soft gluon technique 
has been generalized to the a; — > case in Ref. ^ and one finds that the contribution of 
the gluon g„ can be factorized in fji^ without the need of any coUinear approximation [Q. 
There are three factorized contributions: 



1 



the softest gluon g„ is a real emitted one. One can generalize the result ([l2|) also to the small 
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X case and factorize the colour matrix — J/j„_i (qn) ■ It is surprising that also for x — the 
leading IR contribution for the real emission is given by the eikonal current involving only 
the external partons. Actually, as shown in [Q, for x — the contribution of the emission 
of g„ from the softest exchanged line fc„ is important. This contribution compensate exactly 
the rescaUng k^/(k„ — for the exchanged propagator and, for x — > 0, the emission is 



given by the squared eikonal current ([12D without any coUinear approximation; 



ii) the softest qn is a virtual on-shell gluon not connected to the exchanged gluon k. In this 
case we have again the squared eikonal current factor in (|T^); 

iii) the softest g„ is a virtual on-shell gluon connected to the exchanged gluon k. This factor, 
which is singular only for x -C yn, is positive and proportional to the unit colour matrix. We 

-^e(qf.-k=) (15) 
qn (k + q„)^ ql 

where the last expression is obtained upon azimuthal integration. 

It is an important consequence of coherence in the limit x — that if one takes together 
the first two contributions the colour algebra becomes trivial as in ([l4|) but without any 
coUinear approximation. This can be shown by observing first that in the strongly ordered 
region ( |Tot ) and for x — > one has ?/i — > 1. Thus in this limit the 6 functions in the y variables 
are given by 6{1 — yi) both in the real and in the virtual contribution. Therefore, apart for 
this common 5(1 — yi) function, the sum of the real and the first virtual contribution are 
given by 

fi?„-i(?n)]' {S\k - k«„_, + qj - 6\k - k^„_j} 
^ ^ {8\k - k«_, + qJ - 8\k - k^„„j} , 
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as in the collinear limit. However this expression is vahd in the hmit ?/„ -C Ui without 
requiring any colhnear approximation. To show this observe that for j/„ -C Hi the difference 
of the two delta functions gives a non- vanishing contribution only in the collinear region with 
Oi ~ ynlvfin- Thus for ynlVi — > we can use the approximate expression ( p!^ for the squared 
eikonal current. Therefore it is the coherence of the soft radiation for x — > which allows 
us to compute the eikonal current in the collinear limit. For finite x one can not neglect the 
differences in the arguments of the 6 function in ?/j for the real and the virtual contributions 
and, in this case, the approximation ([l^) is valid only if one requires the collinear limit. 



In conclusion from the factorization results in (15) and ([l6|) we obtain in the soft limit 
Dn y% the recurrence relation in (|TTp which is valid both for x — > (without collinear 
approximation) and for x — > 1 (in the collinear approximation) The three (5-function contri- 
butions in the curly bracket of (|TT|) have the following origin: the first corresponds to the real 
emission; the second and third correspond to the eikonal and non-eikonal virtual correction 
respectively. The new feature of the last contribution, which is relevant only for x — 0, is 
that it involves the total transverse momentum. This implies that the kernel and then the 
solution of the recurrence equation must involve k. This is in agreement with the fact that 



for small x one introduces 111 the structure function at fixed k. 



2.2 General solution of the recurrence relation 

From the recurrence relation (0) one obtains the multi-gluon distributions | f in (|^). 
The virtual corrections exponentiate and can be summed to give (see appendix A) 

-^Vikdv,) |M,(pi,...,p,|^ eL. ^ «^ ^jhov^ K(pi,...,p,)eL., (17) 

1 \ Vi '^Vit 

where the virtual corrections are given by 

Vr{piP2---Pr) = exp<-as — —([[ ^^P]^s ~ ^^^^ 

with 

kj — kj_i p^j , Xj — Xi—i t/i J k — kf , X — Xj. . 

The first exponential in Vr is the eikonal form factor which sums all virtual corrections 
corresponding to the second term in the curly bracket of (pTj). Notice that for y —>■ this 
integral diverges. Similarly all real integrations diverge for ?/j — * 0. These IR singularity 
cancel in inclusive quantities such as the structure function. However, if we want to keep 
separately real and virtual contributions, one should regularize these IR singularities by a 
momentum transverse cutoff. This will be carefully done later. 

The remaining factor in (|T8|) is given by the form factors which sum the contribution 
from the last term in the curly bracket of (p^Ij), corresponding to the non-eikonal virtual 
corrections. These form factors give singular contributions only in the case of soft exchange 
(xj <^ Xj-i), i.e. for x ^ 0. 
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3 Coherence and angular ordering 



The expression of the squared amphtudes | M |^ in (|T^ is given in the strongly y-ordered 
region (p!oD . In order to obtain an evolution equation one has to exchange y-ordering with 
some other ordered variable involving the transverse momentum or angle of emitted gluon. 
The identification of this variable allows us to exhibit the coherence properties of these 
distributions. One can show that within the leading IR accuracy one can replace in (|T^ the 
ordering in y with the ordering in the angles between the incoming gluon and the emitted 
ones, i.e. the ordering in ^ (see Appendix B). We have 

TT— ^ Vr{piP2 ■ ■ ■ Pr) ^1 r =^ IT" T K(PlP2 • • • Pr) ©r...! , (19) 

1 Vi mt 1 Vi mt 

where 

el..i^e(e.-e.-i)...0(6-6), 

so that the gluons {pi . . .pr} are now emitted in the angular ordered phase space region 

6 < . . . < ^r. (20) 

Since the ?/-variables are not any more ordered we introduce the usual Zi variables 

SO that the y-integration factors in (|17[) can be written 

^ f[ ^s{l-x-j:y.) = n dz, as{-^ + -] 6{x - z, ■ ■ ■ Zr) . (21) 

We recognize the singular parts of the gluon splitting function. For ^ 1 we have soft 
emission and for Zi ^ Q soft exchange which is relevant for x — > 0. The finite contribution 
to the gluon splitting function —2 + Zi{l — Zi) can be included. They can not be obtained 
within the soft gluon techniques. 



The IR singularities — > 1 present in the real emission distribution ( PT| ) are cancelled, 
at the inclusive level, by the virtual corrections in V^. If we want to keep separately the real 
and virtual contributions we need to regularize these IR singularity by introducing a cutoff 
Qq in the transverse momenta. As usual we introduce the rescaled transverse momenta 

qi = = ^i-i\fsii, (22) 

i Zi 

so that the angular ordered region (pO]) becomes 

ii > 6-1 , 6 > 6 =^ C[i> Zi^iQi-i , gi > , (23) 



where we have introduced also the coUinear cutoff or jj = ^/s^. The IR singularities in 
(pip for T/j — > (i.e. Zi — > 1) are regularized by limiting the emitted transverse momentum. 
Pit > Qo, giving 

l-z^> Qo/q^ ■ (24) 
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In terms of these variables, the virtual corrections in (|T^) can be written, to leading IR order, 
in the form (see Appendix C) 



Vripm ■ ■ ■ Pr) ^ ^siQ, Zrqr)Y[Asiqi, 2j-igi-l) A{Zi,qi,ki), (25) 

1 

where As is the usual gluon Sudakov form factor in the double logarithmic approximation 



A / N \ l'^^ dq^ f^-Qo/Q ^5 I 
Asiqi, Zi.iq^.i) = exp <; - / — / dz- \ . (26) 

(2i_i(5ri_i) q JO i — z 



Notice that the q integration region (g > Zi^iqi^i) corresponds to the angular ordering 
constraint (^3|). Here the regularization of IR singularity {1 — z) > Qo/q is the same as in 
the real emission phase space. 



The other non-Sudakov form factor is given by 



(27) 



The upper limit in the integration comes from the non-eikonal form factor in (0), while 
the lower limit zqi corresponds to the angular ordering > {p = z Xi^i y/s^ > zqi ) and 
comes from a part of the eikonal form factor not included into (see Appendix C). This 
form factor is different from the Sudakov form factor in various respects. The Sudakov form 
factor is associated to IR singularities for the virtual on-shell gluon thus it regularizes the 
z —>■ 1 singularity in the splitting function of the usual evolution equation. It is a function of 
the local variables at the vertex such as the two angular variables in (pUf). The form factor 
in (^) is associated to a singularity in the soft exchange gluon Zi 0. It depends on the 
energy fraction z^, the angular variable and the exchanged transverse momentum kj. It 
is this k-dependence which gives one of the important feature for coherence in the small 



X region. For x finite ( [27| ) gives a subleading correction which is not present in the usual 
Sudakov form factor even at subleading level |^ . This form factor is related to the Regge 
form factor of the gluon Ajii^Zi, kj) which is found in the BFKL analysis (^. Here the lower 
bound integration over p is bounded for z ^ by a coUinear cutoff fi. Due to this connection 
we shell call the non-Sudakov form factor A{z, q, k) also the hard Regge form factor. 



4 Unified coherent branching 



In conclusion, from (0), ([T9|), (|2lD and ( ^Sf ) we can write the k-structure function in eq. (§) 
in the form (see the kinematical diagram in Fig. 2) 

^(x; k, ko; g, /i) = 5(1 - x) 6\k - ko) As{Q, i2)e{Q -12) + ^. As{Q, z^r) 



X Y[ <j — jdziasP{zi,qi, h) As{qi, Zi^iqi^i) Q{qi - Zi^iqi-i) \ 5{x - Xr) 5^(k - k^ 



i=l 



(28) 
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where Xi = z^Xi^i , kj = — q- + kj_i and ki = | kj |. We have included an incoming transverse 
momentum ko and the colhnear cutoff fi = zqQq which corresponds to the factorization scale 
for colhnear singularities. The splitting function is 

P(z„ g„ k^) = - + ^i£^ . (29) 

1 Zi Zi 

We have included the non-Sudakov form factors A only in the Zi ^ singular contribution 
of the gluon splitting function (plf) since A is regular for finite Zi and gives a non leading 
correction. It is easy to include in P{z, q, k) the finite terms —2 + Zi{l — Zi) which are relevant 
in the region of x not large or small. It is also possible to generalize the equation and the 
splitting function to include the quark contributions. The variable Q has been introduced in 
(I) as the upper bound on the total k (see (||)). To the present accuracy this is equivalent 
to set Q as the upper bound of all emitted transverse momenta. 

The main feature of this branching is angular ordering which is present in the real emission 
phase space (g^ > (see (p!9|)), in the Sudakov form factor (g > ZiQi) (see (p6D), 

and in the non-Sudakov form factor (p > zqi) (see (|27|) ). The fact that the non-Sudakov 
form factor A(zj, qi, ki), relevant only at small x, depends also on the exchanged transverse 
momentum fcj is the reason way we need to introduce the unintegrated structure function. 
The branching distribution in ( pSD and ( P9| ) is the basis for the new Monte Carlo program 
131 which simultaneously takes into account coherence for large x (to double logarithms) 



and for small x (to all loops). 

In the small x region eq. (pSD sum all leading IR singularities, i.e. powers of a^/x ln"~^ x. 
This contributions are obtained when all variable z vanish. Here the lower boundaries for 
the transverse momenta vanishes both in the real emission [q^ > —>■ 0) and in the 

non-Sudakov form factor [p > zqi —>■ 0). We have than that the colhnear singularities 
arising from the vanishing of these lower boundaries generate not only logarithms in the 
transverse momenta, which give powers of InQ/^, but also logarithms in the variables z, 
which give further powers of Inx. Recall that in the present formulations all these logarithms 
of colhnear origin are included without approximations. As we shall discuss later, in the k- 
structure function all colhnear singularities fully cancel so that angular ordering becomes 
equivalent to the multi-Regge region. This is way the expansion of the k-structure function 
contains only powers oi ag/x ln"~^ x and the BFKL equation is obtained. 



5 Unified evolution equation 

The distribution J-'{x, k, kg; Q, /i) in ( p8|) is given in terms of the angular ordered phase space 
and therefore one can deduce an integral or a differential equation in the angular variable. 
See Ref. P, for the case a; — 0. Here we consider the extension to the regions of x both 
small and large. 

For X not small one has that all Zi are finite and, to leading colhnear logarithm, angular 
ordering in (^) is equivalent to qi ordering (g, > Zj_igj_i ~ qi-i)- The evolution equation is 
obtained by differentiating (^) with respect to Q and one obtains the usual AP evolution 
equation with coherence i.e. the leading a5'ln(l — x) contributions are correctly summed 
by taking into account the rescaling factor in the angular evolution variable qi = pa/ (1 — Zi). 
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For X — i> one has that Q can not play the role of an angular variable. Indeed the last 
branching phase space is Q > > ^r-i?r-i so that, by differentiating with respect to Q, we 
get Q = Qr while the upper scale for the angular variable must be ZrQr- Therefore, in order 
to obtain an evolution equation also for x ^ one needs to introduce an additional variable 
q giving this upper scale for the last angle of the emission 



We then introduce the distribution A{x; k, kg; g; Q, ju) defined by ( p8[ ) in which we add the 
constraints q > Znqn- We have then 

A{x] k, ko; g; Q, /i) = 5(x - l)5^(k - ko) A5(g, /i) Q{q - /u) + ^ / 0(g - Zrqr) As{q, Zrqr) 



X n <^ — jdzi asP{zi, qu h) As{qi, Zi^iqi^i) Q{qi - Zi_igi_i) > 5{x - Xr) 5^(k - k^ 



=1 



At g = Q we have 

A{x,k,ko;q = Q;Q,ij) = J^(x, k, ko; Q, /i) 



(30) 



and for g < X /i the integral vanish thus 

A{x, k, ko; g = x/i; Q, jj) = S{1 — x)5^(k — ko) . (31) 
This distribution satisfies the following integral equation 

A{x, k, ko; g; Q,fi) =S{x- 1) 5^(k - ko) As(g, /i) e(g - /i) 

/•O' _ _ cPqdz _ X , (32) 

+ / As(g, zg)e(g - zg) — -— as P{z, g, /c) A{-, k', ko; g; Q, /x) , 

with k' = k + (1 — 2;)q. The evolution equation valid for large and small x is obtained by 
differentiating with respect to the last angular variable 

d 



q'^^A{x,k, ko;g;(5,/i) 

(33) 



dq'^' 



Jx z z [\l — z/+ z j \z z J z 

where k' = k + ^^^^q and the azimuthal integration over the direction of q is understood. 
The z ^ \ singularity in the kernel is regularized by the virtual contribution coming by 



differentiating the Sudakov form factors in (|3^) . Notice that the finite —2 + z{\ — z) should 



be added if one considers regular contribution for z away from the boundaries z = and 
z = 1. This evolution equation can be generalized to include also the quark contributions. 



It is important to notice that ( ^3]) is actually an evolution equation in the angular variable 

q/x = and not in g itself. This is seen in the fact that the derivative of the distribution 
^(x, k, ko; g; Q, /i) in the l.h.s of (|33|) at x and g is given by the distribution at x' = x/z 
and g' = q/z. Moreover the boundary condition ( pi]) is given at the minimum angle y/s^o = 
g/x = /i. Of course one can use instead a boundary condition at fixed g. We have than that 
the solution of (|33|) is obtained by increasing g/x up to the values q < Q. 



12 



Consider the two limits x ^ 1 and x — 0. In the first case we recover easily the usual 
evolution equation for the structure function F{x, Q) in with coherence included This 
is seen by observing that for finite x also z is finite. Recalling that in this case we required 
only the leading collinear accuracy, we can make the following approximations in (pSl): i) 
take as ordering variable q = Q; ii) neglect the non-Sudakov form factor which it is not 
singular for finite z; iii) replace in the integrand q/z with q. With all these simplification we 
obtain 



d dz 
Q'^Tr^^{x,k,ko]Q,fi) = —as 
OLJ J X Z 



/■I dz 




1 


/ —"5 




+ - 


Ix Z 




+ z 



^(-,k',ko;Q,/i) 

z 



with k' = ^+{1 — z)Qn and the integration over the two-dimensional versor n is understood. 
Notice that by neglecting the non-Sudakov form factor, there is no k dependence in the 
kernel. Therefore we can integrate this equation over k and obtain the usual AP equation 
for the structure function (a part for the finite part in z for the splitting function). The k 
dependence of T{x, k, ko; Q, /x) is known [|l| and is of a Sudakov form factor type. 

We consider now ( p3D in the small x region where we can neglect the 1/(1 — 2;)+ term in 
the kernel. We obtain (see also 0, ^) 

<f^A{x, k, ko; g; Q, ^^) = f A(z, ^, k) ^ k', ko; ^; Q, /.) e(Q - ^) , (34) 
oq'^ Jx z z z \z z / z 

with k' = k + — Notice that for Q and xQ < q < Q the 2— integration is over the range 
q/Q < z < 1. Therefore this distribution becomes independent of g as g approaches Q and 
satisfies the BFKL equation. This is analyzed in the next Section. 



6 Coherence at small x 



At small X we can neglect in (|3^) both the real and virtual contributions which are singular 
for z ^ 1. The corresponding approximation in (|30|) gives 

A{x, k, ko; q; Q,fi) = 5{x- 1) 5^(k - ko) Q{q - fi) 

+ XI / ~ ^rqr) n I —^dzi — A{zi, qi, h) Q{qi - Zi-iqi-i) \ 5{x - Xr) 5^(k - k^) , 

r ^0 \ nqf Zi J 

(35) 

with Xi = ZiXi-i , kj = kj„i — and fcj = | k^ |. 

As we shall see, in this fully inclusive distribution all collinear singularities in the real 
emission phase space q^ —>■ Zi^iqi^i — > are compensated by the collinear singularities in 
the non-Sudakov form factors (^) for p —>■ zqi 0. Therefore at the inclusive level one 
can simplify the phase space by substituting angular ordering with the multi-Regge region 
in which all transverse momenta are independent and by regularizing the exclusive collinear 
singularities in ( P3[ ) by the cutoff fi in both the real emission and the virtual integrals. 
After performing the sum the limit /x — is finite. This corresponds to make the following 
approximations in ( ^5|) 

0(g* - Zi^iQt-i) &{qi - ^) l^gg^ 

A{zi,qi,ki) AR{zi,ki) . 



13 



Within this approximation we obtain the following representation of the k-structure function 
^o(a:, k, ko; Q, fi) = 6{1 - x) 6\k - ko) e{Q - fi) + 

J2 [ f[\^^dzi—Aji{zi,ki)\6{x-Xr)6'^(k-kr)' ^^'^'^ 
" =1 -^^^ 1 I ' 



in which the phase space is the multi-Regge region (see 



6.1 Perturbative expansion 

We check explicitly to four loops that the small x structure functions in (|35|) and ( pT)) are 
the same for large Q. We will see that the real and virtual contributions have different type 
of singularities. In the coherent formulations ( |35D we find stronger singularities for x — 
than in the approximate form in (|37[) . All these stronger singularities do cancel in the fully 
inclusive sum (k-structure function), but they do not cancel for associated distributions. 

Consider the contributions to the structure functions for a fixed number r of emitted 
initial state gluons 

F^{Q) = Cdxx^F{x,Q) = l + Y.Fl'\Q). (38) 

In the formulation with coherence the structure function is obtained by integrating (|35|) over 
k and setting q = Q. For simplicity we assume ko = and we obtain 

Fu^iQ, yu) = 0(<5 - Ai) + / n — ^dzi — z'^ A(zj, g^, ki)Q{qi - Zi^iqi^i) , 

^1 Jo Y vrg/ Zi 

with the collinear cutoff yU = z^q^. 

In the formulation without coherence in (^) we obtain 

i^o.(g,/i) = e(g-/i) + 5: / Yl^dz^^z^Aniz^ki). 

^1 J^^2 Y nqf Zi 

For the structure function Fq^{Q) in ( P7| ) without coherence we readily obtain the pertur- 
bative expansion. For large Q one obtains to four loops (T = ln(Q//i)) 

l{2^Tf + U2^Tr - li2^Tr + . . . , 

UJ UJ UJ 



f!^1{q) 


= 2^T 

UJ 


F^liQ) 


= {2^1 

UJ 


FiiiQ) 


= {2^2 

UJ 


F^liQ) 


= {2^2 

UJ 



C(3)r + 



(39) 



The perturbative expansion of these contributions is of the form 

oo — n 
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Summing all exclusive emission contribution one finds the perturbative result of the 
BFKL anomalous dimension up to four loops. The first correction to the one loop anomalous 
dimension is only at the fourth loop. The Riemann zeta functions are obtained from integrals 
of type 

We neglected the collinear cutoff /i since the integral is regular for g ^ 0. Moreover we have 
taken Q as the upper bound both for the emitted and the exchanged transverse momentum 
which is valid for q' <^ Q. 

In the coherent formulation the perturbative contribution to four loops are 

f <..(Q) = ?5£t - <^T + s'-^T - 15 <^T + . . . . 

UJ UJ-^ UJ^ Uj' 

Fi^HQ) = (2^.)^ { . I! + I! + 2-1) - (2^.)* {45 + 6^ + 22^ + 1C(3)^) + ■ ■ • , 

(41) 

The Riemann zeta functions are obtained from the same integrals in (^OD with the lower 
bound /i in general substituted by zq'. The perturbative expansion of F^'\Q) is of the form 



, ,2n—m 
n=r m=l ^ 



The singular terms with m < n are related to coherence since they are coming from angular 
ordering, i.e. the z dependence in the emission phase space {qt > Zi-iqi-i) and in the virtual 
integration {p > zqi). By summing all the exclusive contributions one finds that all these 
more singular terms cancel and obtains the same anomalous dimension dimension (^) to four 
loop. This is an example of the complete cancellations of the collinear singularities in the 
k-structure function. 

At a less inclusive level, such as for the associated distributions, the collinear singular 
terms {a'g/uj'^^'''^) with m < n do not cancel any more. In the approximation (|36D instead. 



the perturbative coefficients are given only by powers (0:5/0;)". Therefore the approximation 
of neglecting angular ordering by simply putting the cutoff /i (see (p6D) is not any more valid 
for the associated distributions. 

As an example we consider the average number of gluons emitted in the initial state, 
i.e. we neglect the final state branching of these gluons. This average is obtained by summing 
the contributions (|3^) and (^l]) with a weight given by r. For the formulation with and 
without coherence we find 



UJ Z \ UJ 
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6.2 Relation with the BFKL equation 



We want to show now that the structure functions at fixed k satisfy the BFLK equation for 
X — i> 0. Consider first the case of no coherence obtained by approximating angular ordering 
with the multi-Regge region. The unintegrated structure function jFo^(k, ko; Q, /i) is given 
in (^). We show that the dependence on Q and is of higher twist type and, neglecting this 
dependence, the unintegrated structure function satisfies the BFKL equation. From ( P7[ ) we 
have 

and one readily deduces the following integral equation 



^oc.(k,ko,Q,/i) = (52(k-ko)e(g-/i)+ / ^rfz^"^A^(^,A:)J-o.(k+q,ko;g,/i), (43) 



as d'q 



uj nq 



2 {J^oc^(k + q, ko; Q,fi)-Q{k- q)J^ouj{k, ko; Q, /j,)} 

1 (44) 



l + ^ln^e(A;o-/i) 



The kernel of this equation is regular for the coUinear limit g ^ so that in the integral 
we do not need the cutoff fi. The cutoff can be removed from the inhomogeneous term by 
defining the unintegrated structure function for which we can take the limit /i/Q — »• 



^oa;(k, ko) = J^o^(k, ko; Q, fi) 



l + ^ln^^e(A:o-/i) 



Q 







(45) 



This distribution satisfies the BFKL equation in (§) with the inhomogeneous term 5^(k — ko). 

We want to show that, due to the cancellation of collinear singularities, the unintegrated 
distribution ^^^(k, ko; g; Q, //) obtained in the present formulation, satisfies to the present 
accuracy the BFKL equation. From (p5|) we obtain (see p|, 01) 



(46) 



^^(k, ko; g; Q, /i) = ^^(k - ko)e(g - /i) 
<Pq 



+ 



dz z"^— A(z, g, k) &{q - zq) A^{k + q, ko; g; Q, /i) 



vrg^ 



Integrating by part 



/' dz z^ — A(z, g, k) e(q -zq) = — {Qiq - g) - ^ dz ^ (Mz, g, k) Q(q - 
Jo z u [ Jo oz 



zq)) 



we obtain (the Q and /i dependence is understood) 
o^s f^^ d'^q 



c^s f d q 

Au;ik, ko, q) = — / — y {Ac^{k + q, ko; q) -Q{k- q) Au;{k, ko; g')} + 5a;(k, g) 



+ 52(k-ko)e(g-/i) 
with q' = min (g, g) and 



1 + ^ ln^e(A;o-/i) 



UJ 



(47) 



Su;ik,q) 



UJ Jq 



as fQ' d^q 



Tcq^ 



A^{k + q, ko;q) 
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A(^,g,k)-1 
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1 -1 



(4J 



As before we remove the n dependence in the inhomogeneous term by defining the distribu- 
tion 

In the equation for this distribution we can remove the cutoff and take the hmit /i/Q 0. 
The contribution 5^(k, g) vanishes for q ^ Q and does not contains leading contributions 
given by the powers a^/oj'^. By neglecting 5(^(k, q) we have that the integral equation for A 
becomes the BFKL equation and in these limits we can approximate 



Aui{k, ko; g; Q, /i) ^ T^{k, ko) 



(49) 



7 Associated parton distributions 



The coherent formulation can be used to compute not only the structure function but also 
the associated distributions for any value of x large or small by generalizing the jet calculus 
IP . As discussed in the previous section the present formulation with coherence is equivalent 
to the BFKL equation only for the fully inclusive case (k-structure function), but different 
for distributions of the associated radiation. 

As illustration we consider the single inclusive parton distribution with a given energy 
fraction y in the hard process of Fig. 1. The relevant distribution is described in Fig. 4 and 
is given in terms of the distribution A{x, k, kg; Q; g, /i) as follows (take ko = for simplicity) 

E{x,y,Q,^i) = [ dx' d'k d'k'^A(^, k, k'; q = Q;Q, zq) / D ( / p,] 
J x'z \x'z J x'[l — z) \x'[l — z) I 

fQ d^q 

X / — -dzasP{z,q,k')A{x',k! + ^^,kQ■,q■,Q,^i), 

(50) 

where = (1 — z)q is the initial state emitted transverse momentum, P{z, q, k) is the 
splitting function in ( ^9]) which includes the non-Sudakov forma factor and D{y/x", q) is the 
inclusive distribution of the final state gluon of energy fraction y emitted in a jet of energy 
fraction a;" = (1 — z)x' and the hard scale is = (1 — 2;)g as required by angular ordering. 
The maximum angular variable is set at g = Q as required by (p2D . From this expression we 
can compute for instance the average associated multiplicity. Notice that in these cases one 
maximises the value of pf 

Consider the region of small x in which all z variables are small. Taking the moments 
one finds 

Si^,a;'(<5,/i) = J dxx'^ dyy'^' T,{x,y;Q, jj,) = J d'^kd'^k' Au,{k,k';q = Q;Q, zq) D^i{q) 

X / -idzz'^{l-zr ^A{z,q,k')A^+^>{k' + q,ko;q;Q,fi), 
where for small x we have approximated pt — q. 

It is easy to check that this distribution satisfies the energy sum rule, namely one has 
S.,i(Q) = jdxx^{l- x) F{x, Q) = F^iQ) - F^+^{Q) . 
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We take u' = 1 in (0) and use the energy sum rule D^i=i{q) = 1. We have 

= /rf'krfVX(k,k';g = g;Q,^g) /"^ ^ dz (z'^ - z^+') ^ A{z, q,k') 
J JiJ? Tiq^ z 

X Aui+i{\i! + q, ko; g; Q, /i) 



(52) 



For the second term we use the small x evolution equation (B^. Changing integration 
variable from g to g = we have 

g'— X+i(k',ko;g;Q,/i) = / dz z'^^^ — A{z,q/ z,\^) A^+i{\^+(\/ zM\q/ z\Q, ^l)Q{Q-q/ z) . 
oq"^ Jo z 

For the first term we use the evolution equation in the minimum angle 

g2— X(k, k' + ci;q;Q,q) = - [ dz ^(k, k'; g; Q, zq) ^ A(z, g, k')e(Q - g) . 

ag^ JO -2 

Finally, changing variables k' + q to k' in this second contribution we have 

^LoAQ) = - J Jq2 '^^^^{■^'^(k)k';g = (5;Q,g) ^aj+i(k',ko;g;(5,/x)} . 

By using the boundary condition 

^^(k, ko; q = n;Q,fx) = 5^(k - ko) , 

and 

F^iQ) = |rf2kX(k,ko;g = g;g,/x), 

we find that the sum rule is satisfied. 

It is interesting to compare (|^ ) with the corresponding expression obtained neglecting 
coherence by extending the multi-Regge phase space [jl6|, |T^. As an example consider the 
case in which we neglect the final state radiation, i.e. we set D^/{q) = 1 in (|5T|). If we take 
uj' = we are computing the average number of emitted initial sate gluons. This average has 
been computed perturbatively up to the fourth loop in the previous Section. This show that 
the high singular terms Og/o;^"'"™' with m < n, present in the formulation with coherence, 
do not cancel. Therefore, in this case the approximation of substituting angular ordering 
with multi-regge region is not valid. 

In general we can identify in (^Tj) the origin of these high singular terms. From (^5]) and 
(^) we have that the coUinear dependence on the cutoff zq in the first distribution A inside 
the integral of (0) is given by a factor 

l + ^ln-^0(A;'-^g). 

id zq 

This factor is not compensated by the integration over the non-Sudakov form factor due to 
angular ordering constraints k > zq and this generates a§/a;^'^~™ with m < n contributions 
not present in the multi-Regge region. The conclusion is that for the associated distributions 
in general angular ordering can not be approximated by multi-Regge regions. 



18 



7.1 Jet calculus at small x 



The expression of the single inclusive distribution in (|5T| ) can be generalized by using the 
jet calculus algorithm. This algorithm is based on the fact that the branching process is 
factorized. For instance the generating functional of y distributions can be obtained as 
follow. Consider first the generating functional Gt[l, q, u] of inclusive distribution in a gluon 
jet of energy fraction 1 and angular variable q which depends on the function u{y). This 
functional is normalized by Gt[l, q,l] = 1 and the single inclusive y distribution is given by 

Diy, Q) = Gt[l, g, u] \u=i ■ 

The general ra-gluon inclusive distributions are obtained by further functional differentiation 
at u{y) = 1. Then we define the corresponding functional for the hard process at scale Q of 
initial state radiation at small x. We define the space like functional k, ko; g; Q, /i; u] 
with the previous meaning for the variables. By using the factorization of the branching 
process this functional satisfies the following evolution equation 

pQ^ (pq ^2 

GJx,k,ko;g;Q,yu;M] = 5(1 -x)5^(k-ko)6(g-/i) + / Q{q-zq) — - — 

Jo 7rg2 z ^53) 

as ] X 

— A(z,g, k) Gt[x',q,U:^>]'> G^-, k + q, ko; g; Q, /i, m] , 



X 



where x' = {l — z)x/z and Ux{y) = u{xy). From p3| ) we have that at u{y) = 1 this functional 
is the unintegrated distribution ^(x, k, ko; g, /x). The single inclusive distribution in (^) is 
obtained by differentiation 

E(a:, y,Q) = [ cfkcfk' ^ Gs[x, k, ko; Q; Q, /i, u] \u=i ■ 
J bu[y) 



To see that this distribution coincides with the one in ( pTl) one proceeds in some formal way 
exploiting factorization which is the bases of jet calculus. We write (|53|) in the formal way 

G,\u] = l+K[u]-Gs[u], 

where the kernel K[u\ corresponds to the distribution in the curly bracket in (pSj). Taking 
the functional derivative one has 

^ Gs[u]\u=i = Gs[u]-^-^-Gs{u]U=i. 



6u{y) 6u{y) 

The expression in ( pOD is obtained from the fact that Gs[l] is into account that Gs[l] is the 
distribution A. 



8 Conclusions 



We recall the main features of the analysis here presented for small and large x. 
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1) At the level of the fully inclusive distributions at small x (k-structure function) the 
dominant region of phase space is the multi-Regge region with strong Xj-ordering and with 
transverse momenta of the same order. This is due to the complete cancellation of the 
collinear singularities which arises when two emitted transverse momenta are very different. 
This property is the basis of the BFKL equation and of the fact that the hard scale associated 
to collinear singularities is lost. 

2) For associated distributions the collinear singularities do not cancel in general and 
the dominant region of phase space is determined by angular ordering. The analysis of 
the associated distributions for small x requires the complete structure of collinear singular 
terms, i.e. no collinear approximations. This level of accuracy can be actually attained for 
X — > by using soft gluon factorization. While for the amplitudes the soft gluon factorization 
does not require any collinear approximation, for the distributions one is usually forced to 
introduce collinear approximation in order to perform the multi-gluon colour algebra. For 
X ^ one avoids this approximation and the basis of this fact is in eqs. (^) and ([T6|) . 
Due to the real-virtual cancellation of soft singularities and to the coherent structure of the 
eikonal current, one finds that the distribution of a soft gluon emitted by a jet of partons 
does not vanish for x ^ only if the jet of partons are confined into an angular cone with 
aperture of order x. Therefore for x ^ we can use collinear estimates. Notice that the 
approximation involved is the soft approximation, not the collinear one. 

3) QCD coherence has the same origin for small and large x and and in both regions 
the emission factorizes. Therefore one can write a unified evolution equation which resums 
all IR singularities which at small x does not involve any collinear approximation while at 
large x is accurate only to leading collinear level. Contributions from non soft radiation 
as well as contributions from quarks can be included in a natural way (see also fl^). In 
this equation, the most important difference between the two regions of small and large x 
is the presence of the non-Sudakov or hard Regge form factor (|27|) for x — 0. For finite 
X this form factor gives a finite correction while for x — this form factor cancels all the 
collinear singularities in the k-structure function. Thus at this fully inclusive level angular 
ordering becomes equivalent to the multi-Regge phase space. For associated distribution 
instead the cancellation does not take place and one can not neglect angular ordering. Thus 
from these distributions one can reveal the coherence of QCD radiation in small x processes. 
To study these associated distributions we have formulated a unified branching algorithm 
for DIS processes which allows one to compute associated distributions in all regions of x. 
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Appendix A 

Take n = m + i and consider the case in which the i softest gluons qm+i, ■ ■ ■ , Im+e are 
virtual while Qm is real 

l/m+e < ■ ■ ■ Vm+l < Vm ■ (54) 

Since the quantities Yr^ and Ilr^ do not involve the £ softest momenta the two corresponding 
(5— functions can be factorized and one can sum over these softest momenta and reconstruct 
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the form factors. From (|TT]) we have in the region 



1 



dy ' d Q I 

1 I Vm+i ^Qm+i ) 



Integrating over the region (p4\) and summing over i one finds 



i=0 

r^Y f I N \ _ fy-^ dy [Q^ dp^ _ [Vm dy fQ^ dp^ ^ 



(55) 



(56) 



where in ( |5^ ) we have Xm = x + ym+i + ■ ■ ■ + 1/m+i ~ Vm- Iterating this procedure to the 
next to soft emitted gluons we obtain the result in ([181). 



Appendix B 

We consider the exchange of y with ^ ordering in (p^). The factor not obviously symmetric 
with respect to yi and is given by 

^X-rX{T{y,,yi+i,ki) , T{yi,yi_i, ki) = exp \ -as f —f -?-) ' (57) 

1 { -^vt y -'k'^ p ) 

where the 0-function corresponds the strong ordering region ([lO|). Here we have used Xi ~ 
yi+i, X = yn+i, and kj = — + kj_i. Since the quantity in { ^7\ j is integrated over we can 
exchange variable names as follows 

r r 

1 perm. i 

where q^^...^^ = q^^ + ■ ■ ■ + q^. 

Consider the case with two real gluons. There are two permutations in (^8)). The first 
with 1/2 -C yi and < ^2 gives yi — I and Xi = x + ?/2 — I/2 so that 

Q2i^l2Tiyi,y2;qi)T{y2,x;qi2) ~ 6^16^2 T(l, xi; gi) T(a;i, x; gi2) . 

The second permutation with the ordering yi <^ 7/2 and ^1 < ^2 gives 2/2 — 1, Xi = x + ?/2 — 1 
and qi <^ q2 so that we can write 

©21021 ^(1/2, Z/i;<?2)T(yi,x;gi2) ~ 6^16^1 T(l, x; 912) ~ 6^16^1 r(l, xi; gi) T(xi, x; ^12) . 

The two contributions have the same form within the leading IR accuracy and we conclude 
that we can substitute y-ordering with ^-ordering 

Ql2T{yi,y2;qi)T{y2,x;qi2) T(l, Xi; gi) T(xi, x; ^12) 
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Consider the case with three emitted gluons r = 3. After singhng out, as in (|58D , the angular 
ordering < ^2 < with 6321 we have six permutations of y ordering configurations. 
Consider the contribution from the fundamental permutation 

©3210123 T{yi, 2/2; qi) T{y2, 1/3; ^12) ^(^3, x; qus) . 

Since ^1 < 6 < 6 and < I/2 < yi, we have yi ~ 1, X2 = x + y3 ~ 7/3 and xi = X2 + 1/2 ^ I/2 
so that, factorizing 9321 6^23, we have 

Tiyi, 2/2; qi) T{y2, y^; qu) T{y3, x; qua) ~ T(l, Xi, qi) T{xi, X2] qu) T{x2, x; qua) ■ 

The same expression is obtained for all other permutations within the leading IR accuracy. 
Consider for instance the contribution 

©321 ©213 T{y2, yi; 92) T{yi, y^; qu) T{y3, x; qus) , 

where ^1 < ^2 < 6 and ?/3 < Z/i < 2/2- We have ?/2 - 1, 2:2 = x + 1/3 ~ 2/3, xi = X2 + yi ^ I 
and gi <^ ^2- After factorizing 6321 ©213; we have 

^(z/2, yi, <?2) T{yi, y^; qu) Tijj^, x\ qua) - T{1, yi, 92) T{yi, X2; qu) T{x2, x; qua) 

~ T(l, X2; qu) T{x2, x; qua) - T{1, Xi, qi) T(xi, X2; qu) T{x2, x; qua) ■ 

For the other permutations one proceeds as before and finds the same expression so that one 
can substitute, to leading IR accuracy, 

©213 Tiy2, Vi, ^2) T{yi, ya; qu) T{ya, x; qua) ©321 ^{1, Xi, qi) T(xi, Xa; qu) T{x2, x; qua) ■ 

The proof for the general case can be done by induction. A general physical argument which 
allows the substitution of ordering with ordering is given in f^. 

Appendix C 

The form in ([T8|) of virtual corrections has been discussed |Q and will be recalled here 
for completeness. From (|TB|) we have 

1 T/ ^ A - dy di ' _ /--I dy fQ' dp^ 

^0+ y h+ i 1 Jx, y Jk^ p^ 

where in the first term we introduced the angular variable ^ of virtual gluon. This first 
term with the singularity for y,^ ^ ^ can be expressed as sum of Sudakov form factors 
which regularize the singularities in the emitted phase space for qi —>■ 0. To this end we 
introduce in the virtual correction the same IR cutoff as in the real emission, namely for 
the virtual transverse momenta we set the cutoff Qq. In the angular ordered region of the 
emitted phase space 

^0 < 6 < ■ • • < 6r < 1 , 
we have qi = Xi-i^/s^i and introduce the virtual momentum p = y^/s^ > Qq. We have 

Jo+ y J(.^ i ^ V y J ^ 

(60) 

Introducing in the first sum the variables y = {1 — z)xi and q = Xiy/s^ and in the second 
the variable y = zxi^i, we deduce the final result in (pSl) . 
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Figure Captions 



Figure 1: Deep inelastic scattering at parton level. The dotted line represents the off 
shell photon q. For the dominat small x contribution the incoming p and the 
outgoing Pi partons are gluons and the recoiling system is a quark- antiquark 
pair of momentum p' . 

Figure 2: Kinematical diagram for parton emission. 

Figure 3: Graphical representation of the elementary vertex factorization. 

Figure 4: Jet calculus diagram for the associated single inclusive distribution of parton q 
with transverse momentum and energy fraction y. The exchanged partons 
/c, k' and k' + q have transverse momenta and energy fractions given by k, a; , 
k', zx' and k' + p^, x' respectively. 
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